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Abstract 



We continue the study of the time synchronization model from larXiv:1 201.2141 . 
There are two types % = 1, 2 of particles on the line R, with iVj particles of type i . 
Each particle of type i moves with constant velocity V{. Moreover, any particle of 
type i = 1,2 jumps to any particle of type j = 1,2 with rates N^ctij. We find 
phase transitions in the clusterization (synchronization) behaviour of this system of 
particles on different time scales t = t(N) relative to N = N\ + N%. 
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1 The Model and Main result 

The simplest formulation of the model, we consider here, is in terms of the particle system. 
On the real line there are Nx particles of type 1 and N 2 particles of type 2, N = Ni + N 2 . 
Each particle of type i = 1, 2 performs two independent movements. First of all, it moves 
with constant speed Vi in the positive direction. We assume further that Vi are constant 
and different, thus we can assume without loss of generality that < V\ < v 2 . The 
degenerate case v i = v 2 is different and will be considered separately. 

Secondly, at any time interval [t, t + dt] each particle of type i independently of the 
others with probability a^dt decides to make a jump to some particle of type j and chooses 
the coordinate of the j-type particle, where to jump, among the particles of type j, with 
probability Here oiij are given nonnegative parameters for i,j = 1,2. Further on, 
unless otherwise stated, we assume that an = a<ii = 0, ai2, &21 > 0. 
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After such instantaneous jump the particle of type % continues the movement with the 
same velocity v^. This defines continuous time Markov chain 

W 2 (0 = (xi\t), . . .,4l(t);x?\t), . . . ,a# a (f)) , (1) 

where x^(t) is the coordinate of fc-th particle of type % at time t. We assume that the 

initial coordinates x% (0) of the particles at time are given. We are interested in the long 
time evolution of this system on various scales with N — > oo, t = t(N) — > oo. 

In different terms, this can be interpreted as the time synchronization problem. In 
general, time synchronization problem can be presented as follows. There are N systems 
(processors, units, persons etc.) There is an absolute (physical) time t, but each processor 
j fulfills a homogeneous job in its own proper time tj = Vjt, Vj > 0. Proper time is 
measured by the amount v j of the job, accomplished by the processor for the unit of the 
physical time, if it is disjoint from other processors. However, there is a communication 
between each pair of processors, which should lead to drastic change of their proper times. 
In our case the coordinates x^ (t) can be interpreted as the modified proper times of the 
particles-processors, the nonmodified proper time being x^(0) + V{t. 

There can be many variants of exact formulation of such problem, see [U HJ E]- We will 
call the model considered here the basic model, because there are no restrictions on the 
jump process. Many other problems include such restrictions, for example, only jumps to 
the left are allowed. Due to absence of restrictions, this problem, as we will see below, is a 
"linear problem" in the sense that after scalings it leads to linear equations. In despite of 
this it has nontrivial behaviour, one sees different picture on different time scales. 

There are, however, other interesting interpretations of this model, related to psychol- 
ogy, biology and physics. For example, in social psychology perception of time and life 
tempo strongly depends on the social contacts and intercourse. We will not enter the 
details here. 

We show that the process consists of three consecutive stages: initial desynchronization 
up to the critical scale, critical slow down of desynchronization and final stabilization. 
Introduce the empirical means (mass centres) and the empirical variances 

1 N t Nt 2 

* ( °(*) = ^E4°(*). sm = ^(xf(i)-^)(i)) 

* k=l 1 k=l 

for types 1 and 2 and their means 

»M = E^W(t), «i 2 (t) = A*i(*)-A*2(*), Ri(t) = ESf(t) 

Our first result concerns the asymptotic behavior of the empirical means. 

Theorem 1 For any sequences of triplets (Ni,N2,t) such that irdn(iVi, JVjj) — » oo and 
t = t(Ni, N2) — > 00 the following statements hold 

, vi-v 2 frit) a 12 v 2 + a 2 iv 1 

*12(£J -> : , — > : 

«12 + «21 t «12 + a 2 \ 
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Assume now that N{ = [ciN], where Cj > 0, C\ + c 2 = 1. Results of the next theorem 
cover all domain of asymptotical behavior of t(N) as iV — > oo. 

Theorem 2 There are the following three regions of asymptotic behaviour, uniform in 
t(N) for sufficiently large N: 

• if^jp- -> 0, then Ri(t(N)) ~ hx 2 t{N), 

• ift = t(N) = sN for some s>0, then Ri(t(N)) ~ h(l - e~* 2S )N, 

• «/^^ -> oo, then Ri(t{N)) ~ WV, 

where the constant x 2 > is defined by the formula see below, and 

2a 12 a 2 i (vx - v 2 f 



h 



><2(«12 + «2l) 3 



2 Proof 

Embedded Markov Chain To prove theorems [T] and H] we will use an embedded Markov 
chain. Consider the continuous time Markov chain f;Ni,N 2 {t) = 6vi,iY 2 (^> u )i defined by (pQ), 
and let 

n(u) < t 2 (u) <■■■ <r n = r n (uj) <■■• 

be the random moments of particle jumps. Then {r n+ i — T n }^ =1 are i. i. d. random variables 
exponentially distributed with mean , yN 1 ,N 2 — (Nia%2 + N 2 a 2 \) 1 . 

We introduce a discrete time Markov chain CN u N 2 (n), n = 1, 2, . . . , 

with the same state space ]R Ari+7V2 . The idea is that the asymptotic behavior of the con- 
tinuous time particle system £,Ni,N 2 (t) can De reduced to the asymptotic properties of the 
discrete time chain CNuNsi 71 )- Indeed, by the Law of Large Numbers 

n 

N x a 12 + N 2 a 21 v ' 

In other words, if n is large, the value n^jv^JVa = n (^i Q; i2 + N 2 a 2 \ )~ 1 is asymptotically 
equal to the "physical" time t associated with the continuous time particle system £jvi,jv 2 (£)- 
Similarly to the empirical mean and empirical variance we introduce, for the embedded 
chain, 

* fc=i * fc=i 



Evidently, Xi(n) = xW(r n ) and Di(n) = Sf(r n ). In the sequel we shall deal with their 
expected values 

fii(n) = EXi(n), di(n) = ED^n), 
and shall need also the notation 

l 12 (n) = fjL^n) - fJL 2 {n), r(n) = E (X^n) - X 2 (n)) 2 

Closed equation for empirical means Here we prove Theorem [TJ 
The following lemma can be checked by a straightforward calculation. 

Lemma 3 The functions /i;(n) satisfy to the following closed system 

//i(n + l) = //i(n) + [«i2 2 (n) - m{n)) +v 1 ]j NuN2 + a 12 (v 2 - ^i)7m,jv 2 
// 2 (n + l) = // 2 (n) + [«2i (ni(n) - /^(n)) + u 2 ] Jn u n 2 + a 21 {v\ - V2)1n u n 2 

For l\ 2 the equation is also linear and closed. Namely, 

l 12 (n + 1) = h 2 (n) + [- (a 12 + a 2l ) l 12 (n) + (t>i - v 2 )} 7tv 1i7 v 2 + (a 12 + a 2l ) (v 2 - Vj) 1 2 Ni ,n 2 
= l 12 (n) [1 - 7jv Xl jv a ("12 + «2i)] + 7iVi,iv 2 1 - v 2 ) [1 - 7^,7^ (a 12 + a 21 )\ , 

and thus we get 

Z 12 (n) = Zi 2 (0)iT + (l-iT)fl, (2) 

«12 + «21 

= J*^ H+ ( Ia(0 )-J^!L.aW (3) 

"12 + "21 V "12 + "21 / 

where R = 1 - 7^,^ ("12 + "21)- 

The variables S(n) = ct 2 \[ii(n) + «i 2 /i 2 (n) also satisfy the recurrent equations 

S(n + 1) = S(n) + jn!,n 2 [«2i^i + "12^2] 

Thus 

S(n) = S(0) + n ® 2lVl a ) 2V2 = 5(0) + (n7Ar 1 ,jv 2 ) • ("21^1 + "12^2) (4) 
From ([3]) and (jlj) the statement of Theorem [T] follows. 

Empirical variances We would like to get closed recurrent equations for di(n). 
Lemma 4 The following identity holds 

E(z>!(n+1) J (x[ 1 \t),xf\t),i = l,...,Ni,j = l,...,N 2 yt<T n ) = 



D x (n) + a 12 ~f Nl 



jv 2 



(D a (n) - D 1 (n) + (^(n) - X 2 (n)) 2 ) - |^i W + 



+2 



JVi - 1 



(7jv!,jv 2 (vi - v 2 ) (X^n) - X 2 {n)) + 7^,tv 2 («i - ^) 2 ) 



and a similar identity for E(D 2 (n + 1) | . . .) can fre obtained by a simple exchange of indices 
1 <-> 2. 



Proof of this lemma is a straightforward calculation. Taking expectations in the above 
formulae we see that in the equations for di(n) the term r(n) is also involved. 
Consider the vector w(n) = (di(n), d 2 (n), r(n)) T . We have 



w(n + 1) = Aw(n) + f(n) + g, 



(5) 



where A is a (3 x 3)-matrix, not depending on n, f(n) is a bounded vector function of n, 
g is a constant vector 

A = £ + 7j v 1 ,jv 2 £, (6) 

/ -an/Nx -a 12 /N! -a 12 /Ni \ 



-OL\ 2 OL\ 2 OL\ 2 

B = B ± +B 2 = I a 2 i — a>2i ol 2 \ 

-2(ai2 + a 21 ) 



-a 21 /N 2 
V ~N~i ~N 2 ~ 



-a 21 /N 2 
iVx A^ 2 N 2 J 



-a 21 /N 2 
a 12 a 21 



f(n) = 2-f NuN2 (m - n 2 ) Zi 2 (n) ~f NuN2 , ~t 



N U N 2 





I +7JVi,iv a 



"12 - a 12 /Ni 
"21 - ct 2 i/N 2 
-(a 12 + «2i) 



9 = 1n u N 2 ( V 1 ~ V 2) 2 ( "217^,^2 



(7) 



Note that A,f(n),g depend on Ni,N 2 and on other parameters of the model. Note 
that for sufficiently large N\,N 2 (if «y are fixed) all components of the vector ~^ NljN2 
are positive. Note also that the matrix elements of A are all positive for sufficiently large 
N 1 ,N 2 . Thus Perron-Frobenius theory is applicable. 



Spectral properties of the matrix A It is easy to check that the matrix B\ has three 
distinct eigenvalues Ai = — (a i2 + a 2 i), A 2 = 0, A 3 = — 2(a 12 + a 2i ). We will study 
asymptotic behavior of the model for the case when Ni = CiN, N 2 = c 2 N and N — > oo. 
Denote A := C\U\ 2 + c 2 a 2 i, that is 7^,^ = (NA) 1 . 

For large Ni the matrix B is a small perturbation of the matrix B\ 



B = B \ H _B 2 it 

1 TV 



5i + — 

AT 



/ -«12/Cl 

-«2l/c2 
OL\ 2 , «21 



-ai 2 /ci 

-«2l/c2 
«12 , «21 



Cl 



c 2 



Cl 



c 2 



-«i 2 /ci \ 
-« 2 i/c 2 

«12 «21 
Cl C 2 
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We will use perturbation theory to get eigenvalues of B. For \\{N) and As(A^) it is sufficient 
to write 

Ai(JV) = -(a 12 + « 21 ) + O (N- 1 ) , A 3 (JV) = -2(« 12 + « 21 ) + O (iV" 1 ) , (8) 
however for A 2 (iV) we will use the result from [2], that 



A 2 (iV) = -(^i? 2 , fc 0) + O^— J, (9) 

where the column vector <p is the right eigenvector of B\ with eigenvalue 0, the row vector 
if)' is the left eigenvector of B\ with eigenvalue 0, and moreover tp'cj) = 1. One can take 

* = ( 1 + — , 1 + — , 1 )/Z, Z= (a 12 + a 21 ) (— + —)■ 

"12 «21 / V"12 «2l/ 




Substituting these values to (Q, we get 

7v — \ \ • 



^(N) = -^ + o(-L), (io) 



x 2 = 2Z- pi + ^l (11) 



where 

x 2 = 2Z- 1 ( 

. Ci c 2 

Denote a±(N), a±(N), <Ti(N) the eigenvalues of the matrix A. From ([6]) and (IE])-© we 
have the following assertion. 

Lemma 5 The eigenvalues of A are 

a 2 (N) = l-h- + o(4 R ) (12) 



iv 2 - \n\ 

for some positive constants b±, b 2 , 63. 

We will need also the eigenvectors of A, which we denote correspondingly by e[ N \ , ejj . 
It is clear that they are also the eigenvectors of the matrix B\ + hB 2 ^. Using the perturba- 
tion theory p] we conclude that e[ , e 2 , are small perturbations of the eigenvectors 
ei, e 2 , e3 of the matrix B\. Thus, calculating e±, e 2 , e3, we get 

It is clear that with (jSJ) and (fTOl) it is not difficult to find explicitly the constants 6j. 
We will need only 6 2 : 

A Z ci«i 2 + c 2 a 2i 
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Some lemmas on the asymptotic behaviour The solution of the equation fl5]) can 
be uniquely written as 

n 

w(n) = A n w(0) + ^"7(" - j) + (1 - A)-\l - A n )g. (14) 
i=i 

The following result shows that the first and last terns in (IT4l) do not influence the asymp- 
totics of w(n). 

Lemma 6 The following estimates hold uniformy in n and N 

\\A n w{0)\\ < Const, ||(1 - A)-\l - A n )g\\ < Const. 

3 

Proof. Using the basis of eigenvectors of A we can write w(0) = kw,i e i • Then 

i=i 



A n w(0) = k wl (a,(iV)) n e 4 (Ar) . Note that sup e[ N) < oo. Moreover, \cn(N)\ < 1, start- 

^— ' ' N 

i=l 

ing from some N. Then the first estimate follows. To get the second estimate we write 

3 



9 = i 2 Nl , N2 Y, k U N y {N ' 



with the coefficients k° gi (N), bounded in N. Apply the operator (1 — A) 1 (1 — A n ) to the 
latter expansion and note that by Lemma [5] 

2 i - (o-mr , const ._ 

1 - a t {N) ~ N ' * 

and 

, 1 - MN)f 



7 ^ l-a 2 (N) ~ C ° nSt 

Then we get the estimate. ■ 

Now we will analyze the second term in ( TT4l) 

n 

i=i 

Note that the vector function f(n), defined by the formula (J7|), is known explicitely with 
the formula (j2J). 

3 

Let £i, £ 2 , £3 £ K. are such that (0, 0, 1) T = ^e^. Then we have immediately that 

i=i 

a 2 i - CK12 t ai2«2i e 1 

si — 7 ; v? ' S2 — 7 ; , £3 



«12 + «2l) 2 ' («12 + «2l) 2 ' («12 + «2l) 2 ' 
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If £i(N) G K is the coefficient in the expansion "(f^i at 2 = 5^&0^0 e £ > then obviously 

i=i 

6(iV)=e i + 0(iV- 1 ), z = l,2,3. (15) 

We have then V N (n) = Y,U N ) (^1n u n 2 («i - « 2 ) E ^(n-j) (^(JV)) 3 '"^ ef 5 . By ([15j), 

and neglecting O (N" 1 ), we have that for N — > oo the asymptotics of V/v(n) coincides with 
the asymptotics of the sum 

V&(n) := X> (2 TiVl)i v 2 ( Vl - V2 ) f>(n-j) {a t {N)) 3 A ^ ■ (16) 

i=i V j=i J 

Lemma 7 For n = N6(N), where 9(N) — > +oo ; the asymptotics of V^(n) is defined by 
the second term, that is the first and third are small with respect to the second. 

Remind that o~ 2 (N) is the maximal eigenvalue of the positive matrix A. Then this 
corresponds to Perron-Frobenius theory. 

Proof. Consider the formula (j2J). If we assume, that /12(C)) < 0, then from v\ < v 2 it 
follows that each coordinate of f(n) is positive for all n (we use it below). Moreover, 
from (j5J) one can get that there exist constants C\ > C 2 > 0, which do not depend on N 
and n, such that 

< C 2 < {v x - v 2 )l 12 {n) < d Vn, N (17) 
Thus, the coefficient of in the sum (TT6l) is positive and can be estimated from below as 

l-a 2 (N) 

Similarly, for i = 1, 3 the absolute values of the coefficients of e-^ in the sum (TT6j) can 

be estimated from above as 



z — ' 1 — 0" 



■ i(JV) 

Thus, to end the proof of the lemma it is sufficient to compare the asymptotics of the 
following three functions 

l-MAO)^ 1 - (a 2 (N)f eW 1 - MN)) N9{N) 
' l-a 2 (N) ' 1 - a 3 (N) 

and to show that the first and the third are small with respect to the second. It is convenient 
to consider separately two cases: a) 9(N) — > 00, 9(N)/N — > 0, b) 9(N) > cN, and use 
Lemma [5j We omit these details. ■ 
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Asymptotic behavior of expectations of empirical variances Now we are ready to 
study asymptotics of the functions Ri and R 2 and to prove Theorem [2j Remind that the 
intervals between jumps of the embedded chain have exponential distribution with the mean 
lN\,N* — (^A) -1 ) then for large N the connection between discrete time n of the embedded 
chain and absolute time t is n ~ t/lN\,Ni — {NA)t. Thus we can take di( (NA)t) instead 
of Ri(t). Remind also that di and d 2 are the first and second components of the vector w 
correspondingly 

Now we can use the lemma which shows that, as t(N) — > 00, the asymptotics of 
w( (NA)t(N) ) coincides with the asymptotics of the vector 

/ (NA)t(N) \ 

6 \2 lNl , N2 (v l -v 2 ) h2((NA)t(N)-j)(a 2 (N)y- 1 \ e 2 N) . 

As ef } = (l,l,0) T + O(A^~ 1 ), then 

(NA)t(N) 

d i ((NA)t)~2S 2 <y NuNa (v 1 -vz) fl h2((NA)t(N)-j)(a 2 (N)) j - 1 . 

3=1 

To find the asymptotics of this expression, we use the representation (J3]), which gives 
l 12 (n) = C' 1>N + C' 2>N R n , C' hN ^ Vl ' V2 , C' 2jN ^C' 2t00 (A ->• 00). 

«12 + «21 

Note that the following estimate holds uniformly in A and n 

00 1 
< (iVA)- 1 \C' 2tN \ Y,R k = (NA)- 1 \C' 2 , N \ ■ — 



1n u n 2 *£c' 2jN R n -i WN))'- 1 



fc=i 

< Const, 



1-R 



since 1 — R = ■jN 1 ,N 2 ( a i2 + «2i)- Consider now the asymptotics of the following expression 

(NA)t(N) 

2^ Nl , N2 ( Vl -v 2 ) Vl ~ V2 (a 2 (N)y- 1 =2t 

< • n/-.r. -4- flln, 



(N A)t(N) s2 -, / I t\t\ \ (N A)t(N) 

- V1-V2 ( W) i-i =26 j_.(!iz^ .i^Mm 

«i2 + a 2 / A ;; A A a l2 + a 2X 1 - a 2 (N) 



By ([12]) and §T5§ cr 2 (N) = 1 - + O (N~ 3 ) and thus the problem is reduced to the 

N 2 

study of asymptotics of 

\ 2 / / I / A \ \ (NA)t(N) s 



2 6 (V1-V2Y N L L (VA) V 

x 2 (a 12 + a 2 i) I V N 2 ) 



Now the theorem easily follows. 
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